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1. Introduction 

QCD in Coulomb gauge, being best suited to examine the Gribov-Zwanziger (GZ) confine- 
ment ideas [1,2], has been the subject of intense research in the last few years. In a series of papers 
[3, 4, 5, 6, 7, 8], which we will briefly summarize here, we have analyzed the behaviour of all 
relevant two-point functions on the lattice and compared them with the corresponding predictions 
of Hamiltonian variational calculations [9, 10, 11], concentrating on the features relevant for the 
GZ scenario. 

As Gribov in his seminal paper noticed, for non-Abelian theories most gauge conditions admit 
several solutions and the corresponding Faddeev-Popov (FP) mechanism is not sufficient to define 
the functional integral beyond perturbation theory. The field-configuration space must therefore be 
restricted to a domain, continuously connected to the origin, where the gauge condition at hand 
always possesses unique solutions. He then showed how, as soon as the fields cross the boundary 
of such region, the ghost dressing function can acquire a singularity; the "no-pole" condition for 
the FP-ghost at non- vanishing momentum is then necessary to implement the restriction. In partic- 
ular, in Coulomb gauge, he argued how such restriction can imply a diverging gluon self-energy, 
motivating its disappearance from the physical spectrum. 

Many issues remain of course in the above description open. Gribov based his conjectures 
on more or less heuristic arguments. Zwanziger later tried to put the whole set-up on a more 
solid basis, while variational calculations, which are viable in Coulomb gauge since they by-pass 
the explicit construction of the gauge invariant Hilbert space [12], did provide some insight on 
the relation of the GZ-mechanism to the Hamiltonian formulation. In both cases, however, ap- 
proximations need to be made; although many authors tackled the problems during the years 
[13, 14, 15, 16, 17, 18], a satisfactory non-perturbative cross-check from lattice calculation was 
hindered for a long time by the presence of strong discretization effects. We have shown [3, 7, 8] 
how for each propagator a mixture of improved actions and separate treatment of their energy de- 
pendence can quite effectively solve such problems, allowing an explicit check of the GZ-scenario. 

As first suggested in [3], the size of discretization effects can be investigated on anisotropic 
lattices, where the time and space like cut-off a t , a s are kept different. In Fig. 1 we show the effect 
of taking the limit a t — > 0, which controls the lattice Hamiltonian limit, on the SU(2) Coulomb 
gauge functional calculated at fixed space-like cut-off, i.e. RG-point. Perturbative determinations 
of the latter [19, 20] are not precise enough and we have re-checked the only non-perturbative 
calculation found in the literature [21]; details can be found in [8], as well as the details of the 
gauge fixing algorithm, which adapts those introduced in [22, 23]. Following the ideas in [3], a 
first anisotropic analysis in SU(3) had been attempted in [24]. 

From the continuum analysis and from our results in [3, 4] we know that in the pure gauge 
sector the static gluon propagator, the static Coulomb potential and the ghost form factor should 
obey: 
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Figure 1: (a): dependence of the gauge fixing functional fjt (a s ) on the anisotropy | = a s /a t at fixed spatial 
cut-off a s . (b): Deviation of the gauge fixing functional from the Hamiltonian limit, Ft (a s ) —Foc(a s ), as a 
function of the temporal lattice spacing a t , together with its leading power corrections. 



where the Gribov mass M ~ 1 GeV and for the gluon self-energy (Oa = D 1 (p) holds. The quark 
propagator, the fermion self energy and the running mass M(|/?|) take the form [7]: 

urn)- — ,"''?> n2 .-, + ""'7' N> , (1 " 2) 

where Z is the field renormalization function, a the energy renormalization function, nib the bare 
quark mass, m x (mb) the chiral mass and m r (rrib) the renormalized running mass [7]. In the follow- 
ing we shall verify such behaviour and determine the relevant parameters. 



2. Results 



2.1 Ghost form factor 

A careful analysis of the ghost form factor in the Hamiltonian limit a t — > shows that its 
UV behaviour agrees with Eq. (1.1), with y g h = 1/2, confirming continuum predictions, and m = 
0.21(1) GeV, see Fig. 1 (a). In the IR going to higher anisotropics increases the exponent fc g h, as 
shown in Fig. 1 (b), where we plot \p\ K "'d(p), with K m the IR exponent for E, = 1, as a function 
of the anisotropy. The limit E, — > °° gives K g h > 0.5, confirming the GZ-scenario. This however 
disagrees with some continuum predictions fc g h = 1 , deriving from the assumption of the finiteness 
of the static ghost-gluon vertex. Whether this is indeed correct and algorithmic improvements 
could change the lattice result is still a matter of investigation. 
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2.2 Coulomb potential 

In Fig. 3 (a) we show |p| 4 Vc(|p|) as obtained from different anisotropics. Fitting the results 
to Eq. (1.1) we get, in the Hamiltonian limit B, — > oo a c = 2.2(2) a, as expected from Zwanziger's 
predictions [25]. 




Figure 3: (a): Infrared behavior of \p\ Vc(p) / (&7to) for different anisotropies (b): Quark field renor- 
malization function Z(|/?|). 



2.3 Quark propagator 

Our calculations were all made on a set of configurations generated by the MILC collaboration 
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[26], see [7] for details. The use of improved actions is crucial to establish the scaling properties 
of the Coulomb gauge quark propagators. This is very similar to the situation in Landau gauge, see 
e.g. [27, 28, 29], whose techniques we have adapted to our case. 

Fig. 3 (b) shows the scaling of the renormalization function Z(\p\ ) for configurations calculated 
at similar bare quark mass, while the RG-invariant functions <x(|j5|) and M(|j?|) are given in Fig. 4. 
Their behaviour agrees with theoretical expectations, see Eq. (1.2). 



28 J x96, a = 0.086 fm, m = 27.1 MeV, IPG 
20 3 x64, a = 0.1 21 fm, m = 31 .5 MeV, IPG 
20 3 x64, a = 0.1 20 fm, m = 47.3 MeV, IPG 
20 3 x64, a = 0.119fm, m = 63.1 MeV, IPG 
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Figure 4: (a): Energy renormalization function a(|p|). (b): Running mass M(|p|). 



Our most interesting results are given in Fig. 5. Analogously to the gluon self-energy ©a(|/?|), 
the quark self energy (Of {\p\) has a turn-over at \p\ ~ 1 GeV, clearly departing from the behaviour 
of a free particle, and diverging in the IR, see Fig. 5 (a); although awaiting confirmation on larger 
lattices, this would in principle extend the Gribov argument for its disappearance from the phys- 
ical spectrum to full QCD. Moreover, as Fig. 5 (b) shows, the running mass M(|/?|) we obtain is 
quantitatively compatible with our phenomenological expectations from chiral symmetry breaking. 
Fitting it to Eq. (1.2) we obtain b = 2.9(1), y = 0.84(2), A = 1.22(6) GeV, m x (0) = 0.31(1) GeV, 
with^ 2 /d.o.f.= 1.06. 



3. Conclusions 

We have shown that the GZ confinement scenario is realized in Coulomb gauge. The ghost 
form factor d(|/?|) is IR divergent with an exponent fc g h > 0.5, which implies Gribov's no-pole con- 
dition and a dual-superconducting scenario [30] ; the gluon propagator satisfies the Gribov formula, 
implying an IR diverging self-energy, and the Coulomb string tension is roughly twice the physical 
string tension. Moreover from the quark propagator we can easily extract the quark self energy 
(Op(\p\), which is also compatible with an IR divergent behaviour, and the running mass M(|/?|), 
which gives a constituent quark mass of m x (0) = 0.31(1) GeV. 

This is in contrast to Landau gauge, where BRST symmetry seems to be non-perturbatively 
broken, violating the Kugo-Ojima confinement scenario [31], while the GZ confinement scenario 
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Figure 5: (a): Quark self energy 0)f(|p|). (b): Running mass M(|/?|) in the chiral limit m\ — > 0; see Eq. (1.2). 

cannot be realized without the explicit introduction of an horizon function, see e.g. [32] for a recent 
review; its physical implications and how these can be related to the presence of dim-2 condensates 
[33, 34, 35] are an interesting issue still debated in the literature [36]. 
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